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Bose— Einstein condensation at reheating 
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We discuss the possibility that a perturbative reheating stage after inflation produces a scalar par- 
ticle gas in a Bose condensate state, emphasizing the possible cosmological role of this phenomenon 
for symmetry restoration. 

PACS number: 98.80.Cq 



I. INTRODUCTION 

At the end of inflation, the Universe was in a very 
cold, low entropy phase. In particular, its energy density 
was still dominated by the vacuum energy of the fields 
which were previously driving inflation. The process of 
converting this energy into a hot thermal bath of matter 
and radiation is known as reheating, and it represents one 
of the most active areas of research in inflationary cos- 
mology. In the last ten years, following the works [|l],§, 
it has been shown that the first stages of reheating can 
be characterized by interesting nonperturbative phenom- 
ena, according to the particular theory one is considering. 
However, all these different possibilities share the com- 
mon feature that they are followed by a more standard 
stage of perturbative particle production leading to the 
formation of the thermal bath. In the present work we 
focus on the possibility that after inflation some light 
boson field x may experience a Bosc-Einstein condensa- 
tion, after production of the corresponding quanta via 
a perturbative reheating stage. This situation may be 
naturally realized if these quanta carry a charge which 
is either exactly conserved or weakly violated by inter- 
action processes with all species produced at reheating. 
Condensation may then occur if a sufficiently large value 
for this charge is released during reheating . 

The first condition is eventually a statement on the la- 
grangian density governing the x fi^ld dynamics. As an 
example we consider a simple model of a complex scalar 
field X whose interaction term is of the form of a quar- 
tic self-interaction monomial, invariant under a global 
U{1) symmetry. Despite of its simplicity we look at this 
model as a prototype of much more complicated and real- 
istic theories, as for example supersymmetric extensions 
of the electroweak standard model, which admit scalars 
as elementary excitations and conserved, or weakly vio- 
lated, charges. 

Specific mechanisms for production of a non zero and 
sufficiently large value for a definite charge at reheat- 
ing have been widely discussed in the literature, aimed 



to drive, for example, efficient ways to produce baryon 
and lepton numbers. Among these, the Affleck-Dine sce- 
nario 1^ has been shown to be able to produce very large 
baryon and lepton asymmetries. In this case scalar par- 
ticles can naturally develop a Bose condensate. 

In view of these considerations, we think that it is an 
interesting issue to study under which conditions a Bose 
condensate may actually form. Section II and III are 
devoted to this subject. In particular, in section II we 
report a fully general analysis of conditions necessary for 
the formation of the condensate, while in section III we 
discuss its actual formation in a cosmological context. 
We finally consider in section IV the possibility that the 
condensate may have an interesting cosmological implica- 
tion, and leads to the phenomenon of symmetry restora- 
tion. This may reintroduce a monopole problem, since 
when symmetry eventually breaks again, once the expan- 
sion of the Universe dilutes the condensate, topological 
defects can be formed. 



II. CONDITIONS FOR A CONDENSATE 

In this section we discuss the conditions leading to 
the formation of a condensate. We consider a distribu- 
tion of quanta of a light complex scalar field x -i whose 
mass is assumed to be dominated by radiative ef- 
fects. We also assume that the initial distribution has 
a particle/antiparticle asymmetry, Q = — n^, and 
that the charge per comoving volume is conserved (or 
only very weakly violated) by the self-interactions among 
these quanta. ^ Without loss of generality, we consider 
Q > . The production of this charge will be discussed 
in the next section. 



*We neglect here interactions with other fields; see the next 
section for a more general discussion of the system in a cos- 
mological setup. 
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The evolution of the system is of course strictly depen- 
dent on these self-interactions, since they determine the 
thermalization time-scale as well as the radiative mass. 
As an example, let us consider the interaction term 



^int — — h 



XI 



(1) 



At lowest order in h, this term is responsible for 2 2 
scatterings. Combining more vertices, processes which 
change the total number of particles n — n^+n^ are pos- 
sible, preserving however the charge Q. The efficiency of 
these processes in establishing a kinetic equilibrium, i.e. 
the thermalization time-scale, will be discussed in the 
next section. Here we consider the asymptotic equilib- 
rium state in presence of a fast thermalization. 

We start studying the case where a condensate does not 
form. The equilibrium distribution functions are there- 
fore given by 



fBE 



(P) 



(2) 



In this expression, E ; 



m: 



, (T) denotes the energy 



of the quanta, while /3 = 1/T is the inverse temperature 0. 
The presence of the chemical potential /x is related to a 
nonvanishing charge Q . \iQ and /i are both positive, the 
minus sign in Eq. (^) refers to particles x , and the plus 
sign to antiparticles x ■ The equality /i^ = — /i-^ = /i 
holds in presence of effective number violating interac- 
tions such as XXXX*^ XX- 

While the distribution function is always regular, 
we must impose /i < to ensure the regularity of . 
The values of /? and can be determined by the two 
covariantly conserved quantities of the system, that is 
the total energy p and charge Q densities 



(2^ 
(2^ 



[/r(p)+/i"(p)] 



{ir (p) - !r (p 



Let us consider the ratio 



R{p/T,mjT) 



Q_ 



(3) 
(4) 

(5) 



Notice that since we deal with radiation this ratio is un- 
affected by the universe expansion, since both quantities 
Q and p'^/'* scale as the inverse third power of the scale 
factor, so that their ratio R is constant. It is easy to see 
that R monotonically increases with p/T. In particular, 
it ranges between and the critical value Rcr (™x/-^) 
= R {m^/T, m^/T) as p ranges between and . We 



^We use natural units c = Ti = ks = 1- 



have thus a simple criterion to understand whether a con- 
densate forms. From the initial conditions for the sys- 
tem at reheating we calculate the quantity R — Q /p'^^'* , 
which is conserved by the interactions among the quanta 
X and X- If this value is below Rcr {m^/T), the equi- 
librium distributions are just of the form yj^, that is 
the condensate docs not form. When instead the initial 
conditions are such that R is greater than Rcr (m^/T), 
Eqs. (|^) cannot represent the final distributions, but an 
additional singular term is present in /^(p) 



fx (P) = / 



BE , 



(6) 



This shows the appearance of a condensate at zero mo- 
mentum. The physical interpretation of this is very sim- 
ple. Loosely speaking, R > Rcr (m^/T) corresponds to a 
case in which a too large charge Q is present to be stored 
in a regular distribution function with the given energy 
density p . The exceeding part of this charge, denoted by 
Qc, gives rise to the condensate. 

The value of Rcr (m^/T) can be obtained from expres- 
sions (||) and (^) in the limit p — > for the distribution 

. The calculation can be performed numerically, or an- 
alytically when m-^/T is much smaller than one. In this 
regime we find 



p = T^ 




(7) 
(8) 



such that the scaling Rcr {rriy^/T) ~ 0.456 m^/T is ex- 
pected for m^/T ^ 1 . We plot in Figure |^ both the ex- 
act numerical behavior and the analytical approximation 
for Rcr (m^/T). As we see, the analytical approximation 
is very good even for as large as 0.1 0.2 T . 




The value of Rcr (m^/T) versus m^/T both nu- 
merically (solid line) and with an analytic approximation valid 
at small arguments (dotted line). Initial conditions such that 
R > Rcr {rrix/T) indicate the formation of a condensate. 
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If the field x is very light, its mass will be dominated by 
radiative effects. In particular, from the interaction (|^) 
among particles in the thermal distributions (||) we can 
perturbatively evaluate the thermal mass for h <^ 1. 
From one loop correction to x propagator one easily find 



= Ah 



and hence 



d3p 



{2ttY2E 



{ir (p) + fr (p)) 



(T) / h 
T V 3 



(9) 



(10) 



Therefore we see that a condensate forms whenever the 
initial conditions result in 



R > 0.26 /i^/^ 



(11) 



In the next sections we will discuss scalar condensation 
in a cosmological context. It is important to stress again 
that the expansion of the Universe does not modify the 
above condition ( pi] ) , at least until the temperature drops 
below the value of the (non-thermal) mass of x- 

We now consider the case R > Rcr (w^/T), so that a 
condensate forms. In general with the distribution func- 
tion for the scalars x given by Eq. (^), both charge Q 
and energy density p are shared by the regular part f^^^ 
and the condensate. The analysis greatly simplifies in 
the limit where most of the the charge Q is stored in 
the condensate. It is worth noticing that this is also the 
case we are more interested in. When a large charge is 
initially produced, the formation of the Bose condensate 
may lead to cosmologically relevant effects. In this case. 



we have 



Qc ^ and Qc ^ Qth, where we have 



defined with Qc and Qth = Q — Qc the charge accomo- 
dated in the condensate and in the thermal distribution 
respectively. Whenever number density is dominated by 
the condensate, evaluation of the thermal mass using Eq. 
(|) gives 



{2nfQc5^ (p)) -2/1 



Qc 



and therefore 



{Qc) ^ (2hQcf'^ . 



(12) 



(13) 



1/3 

This result of course holds when Qc > T. More gen- 

1 /3 

erally, for smaller Qc /T it is easy to check that a very 
good approximation for the thermal mass is given by 



2hQc + -hT^ 
3 



1/3 



(14) 



The total number density can be estimated as follows 



n Qc + nth{Qc,T) 



(15) 



where the thermal contribution nth{QciT) can be nu- 
merically evaluated using the expression /^|' (p) with 
/ij^ = m^, and given by Eq. (p^. In particular when 
the temperature is larger than the effective mass, though 



it is still smaller than Q 



1/3 



{2hQcY'^ <T <Qll^ , 



we have 



nth 



2C(3) 



(16) 



(17) 



Similarly the energy density can be written as 

p ^ m^Qc + Pth{Qc,T) . (18) 

where pthiQc, T) and m^Qc are the thermal contribution 
and the energy density stored in the condensate respec- 
tively. 
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FIG. 2. The ratios Qth/Qc (solid line), rith/Qc (dashed 
line), pth/m^Qc (dotted line) and R (dash-dotted line), in 
logarithmic scale, versus Qc/T'^ ■ Qth is the charge stored in 
the thermal distribution. Results are shown for h = 0.01. 

In Figure 2 we plot the values for the ratios Qth/Qc^ 
nth/Qc, Pth/rn^Qc, and R, as functions of the variable 
QcjT^i where Qth is the charge stored in the thermal 
component. Results corresponds \,o h — 0.01 and have 
been evaluated using Eq. (|lj). Notice that for Qc ~ 
the energy density is still dominate by the thermal 
distribution, while both the charge and number density 
receive their main contribution by the condensate. 

Changing the value of h slightly affects all quantities, 
but it does not change the monotonic decreasing be- 
haviour as function of Qc/T^ shown in Fig. 2. As an 
illustration we show in Fig. 3 how the above quantities 
depends on h, at Qc/T^ ~ 1, choosing h in the range 
10-10-10-1. 

It is rather important to notice that the total energy 
density p always redshifts as the one of radiation. While 
this is obvious for the thermal part of the distribution, 
for the condensate it follows from the dependence of 



3 



on Qc and T given in Eq. (p^. As a consequence, the 
energy density stored in the condensate evolves in time as 
Qt^^ , where a is the scale factor of the Universe. 
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FIG. 3. The ratios Qth/Qc (solid line), rith/Qc (dashed 
line), pth/m^Qc (dotted line) and 7? (dash-dotted line), in 
logarithmic scale, versus log h, at hxed Qc/T^ = 1- 

As we have said, we will be mainly interested in the 
following in studying the case where Qc represents the 
dominant contribution to n. This condition is satisfied 
whenever Qc > nth- Using the behaviour of nth and p 
with Qc/T^ this condition can be cast in a lower bound 
for the ratio Q/p^/" 

^>--- (19) 

as h varies in the range 10^^" 10^^. Condition ( |l9| ) is 
of course more severe than Eq. ([l^) , which merely states 
the formation of a condensate regardless of its magnitude 
with respect to nth- The order of magnitude of this re- 
sult, obtained by a numerical study, can be also achieved 
analitically by noticing that for all values of the coupling 
we have considered, when nth ^ Qc the energy density 
receives the dominant contribution by the thermal part 
of the X — X distribution, so that 

p^Pth = ^T^ - (20) 
15 

Using the fact that p'^l^ inth (see Eq. (p7|)) we obtain 
that Qc dominates the number density when Qc > nth ~ 
p\l^ /i, i.e. whenever R > 1/3. 

III. FORMATION OF THE CONDENSATE 

In this section we discuss a possible scenario for the 
formation of a condensate in the early stages of the Uni- 
verse, n We start with a field (j) which is oscillating 



For related works, see 




around the minimum of its potential V{(j)). ^ One in- 
teresting possibility is that is a modulus field of a su- 
persymmetric theory. SUSY models offer several candi- 
dates for moduli fields, since their scalar potential almost 
unavoidably presents flat directions if supersymmetry is 
preserved. Most of these directions are expected to ac- 
quire a mass of the order of the supersymmetry breaking 
scale. Actually this scale, in the early Universe, can be 
very different from the present one, due to the dynamics 
of the inflaton or of other scalar fields during inflation 
and the first stages of reheating. See for detailed 

discussions. Nevertheless, though supersymmetry natu- 
rally offers a theoretical framework, we aim to discuss 
the condensation phenomenon in a more general contest. 
For this reason we leave the mass as a free parameter. 
Also for the initial amplitude 0o of the field we have some 
natural scales, as for example the Planck mass, emerging 
in the context of supergravity models, or the breaking 
scale of GUT theories. However, also for (/)o very differ- 
ent values can be considered in principle, and we will also 
take it as an input parameter. 

At very early times the field cf) is frozen at 0o i due to 
the friction provided by the expansion of the Universe. 
When, at a time t = tg , the expansion rate H falls down 
, the field (j) starts oscillating around the minimum 
= 0. The amplitude of these oscillations decreases as 
Q,-3/2 jj: ^Yxis stage is already radiation dominated, we 
then have 

where oq = a (<o) • 

If (p has a sufficiently long lifetime, its coherent oscil- 
lations will dominate over the thermal background for 

> O'dom = ^0 (3/87r) (Afp;/(/)o)^ • After this time, the 
energy density stored in the oscillations redshifts with a 
matter-like behaviour, so that 

/Stt" (j) f a \ "^^"^ 

ff(0 = ^T™*M^(^j ' ""^^'^ ' (22) 

while the amplitude (j) {t) is still given by Eq. (|l|) . 

When eventually the field decays, it produces a re- 
heating stage. In particular, baryon and lepton asymme- 
tries can be generated in a rather natural way through the 
Affleck-Dine mechanism if evolves along a complex 
direction and has baryon and lepton violating couplings. 
The same mechanism can be applied for the production 
of a generic charge Q . Let us indeed suppose that the 
field mainly decays at the time td into the complex 
scalar x considered in the previous section. Through a 



^We assume that the minimum of the potential corresponds 
to = and that anharmonic terms in the Taylor expansion 
of V{(l)) can be neglected, so that has a constant mass . 
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mechanism a la AfHeck-Dine, these decays produce an 
asymmetry 



(td) = P4, (td) /m^ = 0d "^0 



(23) 



Without entering into the details of the production, 
which is cruciaUy dependent on the specific couphng be- 
tween (p and X ^-iid on the dynamics of the scalar field 
(j>, we have parametrized with / in Eq. (|2^) the charge 
produced per quantum of (f) . We have also parametrized 
with the amplitude of the oscillations of (j) at decay. 

It is immediate to see how large / should be in order to 
produce a condensate in the x quanta distribution. Since 
Ql(?l^ — f p,f,{td)^^^ /m^ , we see that the 



condition for a condensate, Eq. (pi]), is met for 



necessary 



/2> 0.07/1^ 

(Pd 



(24) 



The more stringent requirement that the number density 
of X is dominated by the condensate, Eq. (p^), rewrites 
instead 



/' > (0.1 



0.25)^ 

(Pd 



(25) 



The interaction terms responsible for the decay of (j) 
into X will also give a mass to x—x quanta through radia- 
tive corrections. This mass vanishes once 4> has decayed, 
so that the considerations of the previous section (for ex- 
ample, Eq. (p^) remain unchanged. However, it could 
prevent from decaying, whenever {4>d) > 
{we do not consider here nonperturbative decay process, 
where this bound may not apply). This gives an upper 
bound on (jjd , which depends on the specific interaction 
terms between <j) and x ■ For interactions of the form 
cripx^ + h.c, one finds for the decay rate 



r (^ xx) 



ml{(j)) 



16 7rm^ 



^m\-^m\ (0) 



(26) 



which, equated to the expansion rate of the Universe, 
gives 



where we have defined 



0.007- 



'Pi 



(27) 



(28) 



We are assuming that the scalar </) decays after it domi- 
nates. Using Eq. (p^, this translates in the bound 



< 



^0 

\Mpi 



(29) 



We show in Fig. |^ how the ratio (j)d/'m<ji changes for differ- 
ent values of and a . As it is clear from the above dis- 
cussion, see for example Eq. (^), this is a crucial quan- 
tity, since an increase of (j)d over corresponds to an 
increase of the number of produced quanta over their en- 
ergy. We first note that, for any fixed , there exists 
only a finite interval for a where (pd is enhanced. This is 
easily understood, since both a very low and a very high 
a reduces the (p decay rate (^, so that when the decay 
eventually takes place the (p amplitude has been strongly 
damped by the expansion of the Universe. The decrease 
oiT{(p — > xx) with decreasing a is rather obvious, since a 
measures the strength of the x~4' coupling. In the limit of 
high cr, Eq. ( [2^ ) reduces instead to (t>d/m^ ~ m^/ (4ct) . 
This shows that, in this limit, the <p particles are pre- 
vented from decaying by the high effective mass {<p) , 
and only once <pd is sufficiently reduced the decay can 
occur. 




-12 -10 
log (m^/Mpi) 

FIG. 4. Contours for the ratio pd/m^ in the log{m^/Mpi) 
- log{a /m^) plane. Continuous lines correspond, from right 
to left, to pd/m^ = 1, 10, 10^ 10^ and lO". 

For a given , the maximum of pd corresponds to r = 

1/V5, that is for cr, such that ai./m^ ~ A{m^/Mpi)^^'^ . 
From Fig. ^ it is possible to understand the range for 
the coupling between the two species such that a sizeable 
condensate would appear. If a is not too far from this 
value, we have from Eq. (UTt) 



0.05 



Mpi 



1/3 



(30) 



that is (pd/m^ increases as decreases, as also clear 
from Fig. |[ Therefore, formation of a condensate is 
favoured at small tti^ . For a close to the above value, 
the necessary condition Eq. ( p4[ ) , and the more stringent 
condition of Eq. (|25| ) can indeed be cast in the form 



f Z 1-4 



\Mpi 



1/3 



and > (2 ~ 5) 



/ rriA 



PI 



1/3 



(31) 
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However, both the absolute value of (t>d and of Q ~ 
/ (jy^ decrease as becomes smaller. Hence, de- 
pending on which cosmologically scales one is interested 
in, the mass of the field </> cannot be taken arbitrarily 
small, if the condensate is supposed to play any relevant 
role. We will discuss this point in more details in the 
following section. 

Until now we have simply assumed that, when all nec- 
essary conditions are fulfilled, the condensate forms due 
to a fast thermalization of the x quanta. The actual gen- 
eration of the condensate can be discussed by studying 
the typical timescales of self-interactions of the x parti- 
cles. [] 

Kinetic equilibrium among the quanta of x is largely 
governed by 2 ^ 2 scatterings, allowed at tree level by 
the interaction term (|^). At the initial stages, when the 
X are produced by the (j) field the corresponding rate can 
be estimated to be 



(32) 



where n > Q is the number density of the particles in the 
initial state (initially n ^ given in Eq. (p3[)) and (E) 



their average energy, initially of the order of m^. How- 
ever, a complete thermalization requires that also pro- 
cesses which violate the particle numbers are in equilib- 
rium. Among these, the ones which appear at the lowest 
order in the coupling h are 4 2 processes, whose rates 
can be estimated as follows 



Ti_2^h^ny{E)' ,r2 



h'^ n/{E)' 



(33) 



For a thermal distribution without a condensate, one 
has the simple relations n ~ (E')^ ~ . However, we 
will always consider systems for which n > (E)^ . In the 
initial configuration, this occurs for 



<Pd > 



(34) 



Actually, we typically require 4>d ^ , i.e., as we have 
seen, small couplings a < m^, and ^ Mpi, see Eq. 
(p7|). The initial distribution of x , X is thus charac- 
terized by high particle density and small energy of the 
individual quanta, so that the first stages of thermaliza- 
tion proceed via particle fusion. The decrease of n and 
increase of (E) have the general effect of decreasing the 
interaction rates. However, the condition n > (E)^ can 
be never violated if the system has a sufficiently high 
charge Q > p^l^ . This is equivalent to require that most 
of the particle in the final distribution will be stored in 
the condensate, Eq. (p5|). When this condition is met, 
all the interaction rates are thus always greater than the 
ones the x^ X would have if they were in a thermal dis- 
tribution with the same energy ~ ™0</'d ' ^^'^ with- 
out any condensate. If we therefore evaluate all relevant 



rates in the latter case, this will be a sufficient condi- 
tion to ensure that the system we are considering rapidly 
approaches its equilibrium form. 

Among the interactions we have considered, 2^4 
processes have the lowest interaction rate. For a thermal 

distribution with temperature T ~ f 0^ j , we can 

estimate r4^2 ~ /i"* T , so that thermalization is ensured, 
a fortiori, for 



1/4 



M, 



PI 



(35) 



If the coupling h satisfies this bound, the x i X system is 
expected to quickly thermalize right after it is formed at 
the decay time t^, leading to the formation of the con- 
densate, mainly populated by x quanta, and a thermal 
distribution of both x and x- 

For h <^ thermalization can occur only when the 
scale factor of the Universe becomes a ~ a{td) {h^,/h)^, 
since the ratio T/H increases linearly with a. There- 
fore in this case the final effects of the condensate is ex- 
pected to be strongly damped by the dilution of Qc, due 
to the expansion of the Universe. The case h ^ re- 
quires more care, since (as we remarked) sclf-intcractions 
between particles which are forming the condensate are 
much more efficient than in the thermal case {n 3> (E)^) 
and thus the formation of the condensate may be ex- 
pected to hold in this case as well. However, we believe 
that a more precise answer in this regime cannot be ob- 
tained without a numerical study of the relevant Boltz- 
mann equations. For this reason, hereafter the relation 
h>h^ will be always assumed. 

Before concluding this section, it may be convenient to 
summarize the above results. We have considered some 
light scalars x i X with an initial asymmetry Q . Eq. ( [2^ ) 
is the necessary condition for a condensate to appear in 
the distribution of x ■ The distributions for x and x also 
keep a regular part (§), which contributes to the energy 
density of x quanta. The total energy density redshifts 
as the one of radiation, as it is clear from Eqs. (|l8[) and 
(p^). The condition ( p5| ) corresponds to the case where 
most of the particle in the final distribution are stored in 
the condensate, so that the final number density approx- 
imatively coincides with Qc ■ Of course this condition 
is a stronger requirement than just the formation of the 
condensate of Eq. (p^. In the specific example we have 
discussed, the quanta of x are produced by the decay of 
a massive scalar . We parametrized the energy den- 
sity and the decay time of 4> by its classical amplitude 
(f>d when it decays, given by Eq. (p7|). We have assumed 
that (f> decays after it dominates, which is ensured by 
Eq. (^9|). Finally, (^5|) is a sufficient condition to guar- 
antee fast self-interactions among the quanta of x , X j so 
that the condensate actually forms right after the decay 
oi6. 



**For recent studies in thermalization, see fllCMl 
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IV. COSMOLOGICAL IMPLICATIONS 

The aim of this section is to discuss the possible im- 
phcations of the condensate of x particles. In particu- 
lar, we show that it may lead to symmetry restoration. 
As it is well known, this issue has crucial consequences 
for cosmology, since when symmetry eventually breaks 
again, once the expansion of the Universe dilutes the 
condensate, topological defects can be formed, reintro- 
ducing a monopole problem. The possibility of symmetry 
restoration after nonperturbative creation at preheating 
has been emphasized in several works IT^IT^. Here we 
show that this may happen also in the ordinary, per- 
turbative reheating stage, which is generally expected to 
occur in all models of inflation. 

We suppose that the field x is coupled to some other 
scalar field 'ip . For definiteness we will consider -0 to be 
complex, though our result applies to a real scalar field 
as well. Besides this interaction, the field is assumed 
to have the potential 



and to be already settled down in a state with | (V') | = "00- 
Due to the coupling to the i) field, the scalar x ac- 
quires in general an effective mass (?/;) and a modi- 



fied dispersion relation. 



\{Qc.T) 



due to forward-scatterings, which however do not mod- 
ify the distribution function of the x quanta as long as 
m^{Qc) » m^(V'). On the other hand this coupling also 
results in a contribution to the V' one-loop effective po- 
tential MM, 



dAV, 



4, 



1 



d^p6{p^-ml{Q,,T)~ml (0)) 



{fx (P) + h (P)) 



(37) 



In case the x number density is dominated by the con- 
densate, fxip) ^ {2tt)^Qc5^ (p), this relation gives 



\{Qc,T) 



(38) 



where as discussed before we are considering Qc ^ 

This term may shift the ground state of -0 to zero, thus 
restoring any symmetry which may be spontaneously 
broken when instead \{ip) \ = 0o- 

As an example, we consider an interaction term of the 
form g^ |xP- In order to use the analysis contained in 
the previous sections, we first have to constrain our study 
to the case where the x quanta mass is still dominated 
by the radiative term of Eq. (p^, i.e. 



From Eq. (BSl) we therefore get 



Qc 1 



(39) 



(40) 



which shifts the minimum of the total potential 
AV^ to = for 



(41) 



The charge per unit volume Qc takes its largest value at 
the decay of the (j) field introduced in the previous section, 
and is then diluted by the expansion of the Universe. To 
see whether the symmetry broken by ip is restored we 
evaluate the conditions ( p9|) and ( ^ ) at the decay time 
trf. Using the expression (23) for Q^, we have 



00 



(m^ <j)l) 



< min 



(2fe/) 
9 



1/3 



27/6^1/2/^1/6 



(42) 



For definiteness, and to get the largest scale which can be 
affected by the formation of the condensate, we consider 
the case where the ratio ^d/m^ is the largest compatible 
with Eq. (p7|), for each m^. This is equivalent to the 
choice (7* /to0 , so that from Eq. (pO) we have 



(36) 00 



< 



M 



PI 



7770 



2/9 



0.2 



(hf) 



1/3 



0.06 



K1/2/1I/6 



(43) 



Once some specific values for the parameters are cho- 
sen, Eq. (^) indicates which is the maximal scale 0o 
which can be restored by the condensate. It is not un- 
conceivable that scales as high as the GUT one can be 
restored, provided 777^ is taken sufficiently high and the 
coupling a is fixed to give the maximal 0^ as in Eq. (|3^). 
On the contrary, if we are interested in restoring a sym- 
metry at a smaller scale, more freedom is allowed for 
7770 , as well as for the value of (f>d , as it is clear from 
Eq. (^). We remind that these conditions must be sup- 
plemented by the ones found in the previous section, es- 
sentially Eq. (|2j), necessary condition for the condensate, 
and Eq. (^5|) , implying a rapid thermalization of the sys- 
tem. It is worth noticing that these constraints are all 
mutually compatible. 



V. CONCLUSIONS 

In this paper we have considered the possibility that a 
light scalar field, whose quanta are excited at a perturba- 
tive reheating stage, may produce a large Bose-Einstein 
condensate. This phenomenon takes place if these quanta 
carry a conserved charge and the reheating stage pro- 
duces large value for this charge. The motivation for this 
study relies on the fact that in scenarios as the Affleck- 
Dine models for baryon an lepton number generation, the 
decay of a modulus field may typically result in a quite 
large charge per comoving volume in the final relativistic 
degrees of freedom, and this charge is conserved or only 
weakly violated by the self-interaction processes after the 
reheating stage. 



7 



We have studied in details the conditions which should 
be fulfilled for the formation of the condensate for a com- 
plex scalar field with a simple quartic self-interaction 
term, and how the condensation process depends on the 
time scale for thcrmalization via scattering and number 
changing processes. For a wide range for the mass of the 
field whose decays lead to reheating and for the relevant 
coupling constants, the condensation phenomenon may 
easily occur and, moreover, the condensed particles can 
represent the main contribution to the total number per 
comoving volume. 

We have finally considered the role that the condensate 
may have in restoring a symmetry which was broken at 
a larger energy scale. It is interesting to notice that also 
a grandunified symmetry, with scale Mqut ~ 10^^ 
10^^ GeV, can in principle be restored, although this 
requires some fine tuning in the parameters of the model. 
Restoration can occur more naturally for smaller scales, 
until the expansion of the Universe dilutes the charge per 
comoving volume which is stored in the condensate. This 
possibility would represent a new intriguing scenario in 
the evolution of the primordial Universe. 
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